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ALEXANDER WEISS 

Abstract. We propose a simple model for the behaviour of longterm investors 
on a stock market, consisting of three particles, which represent the current 
price of the stock and the opinion of the buyers, respectively sellers, about 
the right trading price. As time evolves, both groups of traders update their 
opinions with respect to the current price. The update speed is controlcd 
by a parameter 7, the price process is described by a geometric Brownian 
motion. We consider the stability of the market in terms of the distance 
between the buyers' and sellers' opinion, and prove that the distance process 
is recurrent/transient in dependence on 7. 



1. Introduction 

In this article we suggest a simple model for the behaviour of longterm investors 
on a share market. We observe the evolution of three particles. One of them 
represents the current price of the share, the second one the opinion of shareholders 
about the share's value, and the last one the opinion of potential buyers. As 
longterm investors do not speculate on fast returns, it is reasonable to assume two 
features: first, the value of the share in the eyes of their holders is much higher than 
the current price, and it is much lower in the eyes of potential buyers. However, 
both groups of investors will not wait forever. They will modify their opinions in 
dependence of the price development. But, as second feature, the traders will only 
slowly adapt to price changes. As opposed to short-time traders, who gamble on 
returns on short time intervals, there is no need for longterm investors to react on 
small fluctuations. 

Eventually, as the price changes and the investors adjust their opinions, the price 
will reach the value, which is expected by the traders. We assume a symmetric 
behaviour of buyers and sellers, and need to consider, what happens if the price 
reaches the right value in shareholder's opinion. Because the price has reached a 
fair level, the investors will sell their shares. At the very moment there are new 
holders, namely the buyers of the shares. Eventually, the price will drop, and there 
is again a group of individuals not willing to follow this decrement. This means, 
while the individuals in the group of longterm investors will change in time, the 
group itself will persist. Figure [1] shows an example for the evolution of the system 
on a logarithmic scale. The price is denoted by B, the opinion of buyers by X, and 
the one of holders by Y. 
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Figure 1. The price B (black), and the opinions X (red), and Y 
(blue) on a logarithmic scale evolving in time. 

We will be interested into the evolution of the distance between X and Y. In 
illiquid markets, i.e. in markets wherein there is only few supply, already smaller 
demands can only be satisfied in connection with a strong change of the price. 
Thus, a large group of traders willing to trade for a certain price provides some 
resistance against further evolution of the price into this direction. Consequently, 
it is of great interest how longterm investors adapt to strong price changes since 
they are providing resistance on levels which are normally on some distance from 
the price. If these investors react to slowly, the price can fluctuate between these 
levels without much resistance, leading to strong volatility. The theory of trading 
strategies on illiquid markets is a very active field of research and there are many 
different approaches to model these markets and their reactions on trading [H El [8] . 
However, the question if large orders on illiquid markets can destabilize them, seems 
to be open. 

Bovier et at describe in [J] a class of Markovian agent-based models for the 
evolution of a share price. Therein they present the idea of a virtual order book, 
which keeps track of the trader's opinions about the value of the share, irrespective 
of whether they have placed an order or not. For practical purposes the model is 
stated in a discrete time setting and in every round one agent updates his opinion. 
As a main feature, the probability to be chosen depends on the distance of the 
agent to the price. In particular, in a market with ./V traders and current price p 
the probability for agent i with current opinion pi to be chosen is given by 



The function h is assumed to be positive and decreasing, reflecting the idea that 
traders with opinions far away from the price react slower to price changes. The 
model is stated in a very general setting, but the authors are able to reproduce on a 
qualitative level several statistical properties of the price process, sometimes called 
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stylized facts, by choosing 
(1.2) h(x) ' 



We pick up on this choice for our model. The logarithmic price process B will 
be a Brownian motion, whereas the opinions of buyers, X, and sellers, Y, are 
described by ordinary differential equations in dependence on parameter 7 > and 
the Brownian motion B. 

The buyers opinion at time t is given by the solution of 

(L3) dt f{t) = (1 + B t - f(t))^ 

whenever X t < B t . By the argumentation above that the individuals within the 
group may change, but the group of traders itself remains, X can hit B, but it is 
not allowed to cross it, and thus, it describes the same movement as B, until B goes 
up so fast that it cannot follow (observe that 1 is an upper bound for the speed 
of X). This happens immediately after the two processes have met, because B is 
fluctuating almost everywhere. As soon as the distance is positive, X is driven by 
(|1.3[) again. Since B is differentiable almost nowhere, some work is needed to give 
a rigorous construction of this process. 

For the opinion of shareholders Y t we assume the same construction with a 
changed sign on the right hand side of (|1.3p . — B is also a Brownian motion, and 
thus we can define equivalently 

(1.4) Y(B) = -X(-B). 

Notice that the speed of adaption to price fluctuations is governed by the pa- 
rameter 7 in our model. Therefore we are interested in the longterm behaviour of 
Y — X as a function of 7. In particular, we would like to know, when Y — X is 
recurrent, and when it is transient. A heuristic argument suggests that 7 = 1 is a 
critical value. For a constant, c > 0, we scale time by c 2 and space by c. We denote 
the scaled versions of the processes by adding superscript c. By Brownian scaling 
we have that B c is equal to B in distribution. On the other hand, X c solves 

(L5) n x * = {i/c+m-xtv 

If one assumes B^ — X^r to be larger than 0, the slope tends to infinity for 7 < 1, and 
to for 7 > 1 as c becomes large. This observation suggests that Y — X remains 
stable for 7 < 1, only. In this paper we show that this first guess is right, and prove 
a rigorous statement about the stability in dependence on 7. 

The remainder of this article is organized as follows. In Section [2] we define the 
particle system formally. X , or Y respectively, will be constructed pathwisely as 
a sequence of processes. The existence of these limits is stated in Lemma T2.ll its 
lengthy proof is given in Appendix [X] In Section [3] we present the main theorem 
and its proof, and in SectionQ]we will discuss, what our results mean for the opinion 
game from 0]. 

2. Construction 

We introduce the processes B, X, and Y formally. While B = {Bt) teR + is just 
a Brownian motion on a probability space {Q, & , (=^t) teK + , P}, X is constructed 
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pathwisely by introducing a sequence of random step functions B e (uj) 1 for which 
the distance to B(w) is uniformly smaller or equal than e. The construction of X e , 
attracted to B e in the sense as explained in the introduction, turns out to be easy. 
At last, we show that X e has a limit as e tends to zero, and call this limit process 
X. The construction of Y follows immediately afterwards. The advantage of a step 
function approach is the simple transition to a discrete setting, which we will use 
extensively in the proof of the main theorem later on. 
For any e > wc define jump times by <7q := and 

(2.1) cr| := min jt > a^ x : B t - B sti > e| ; i £ N, 

neglecting the e-index, provided no confusion is caused. Furthermore, we define 
step functions B e : [0, oo) — > R by 

(2.2) Bt :=B St fort £ [a u a l+1 ) . 
Observe, by definition 

(2.3) sup | S t - Bl\ = e a.s., 

t>o 

and thus B e converges to B on the whole [0, oo) in sup-norm. As already mentioned 
in the introduction, we basically want X to fulfil 

(2.4) ±X t = {l + B t -X t )-\ 

as long as X t < B t . If we substitute B by a fixed number b > 0, the ode (|2.4p is 
explicitly solvable. The solution of 

(2-5) j t f(t) = (l + b-f(t))->; /(0) = 

is 

(2.6) h(t, b) := b + 1 - ((6 + 1) 7+1 - ( 7 + 1) t) ~ . 

We will call h(t, b) well-defined if 

(b+ - 1 



(2.7) b > and t < 



7 + 1 



Obeserve that the bound on t ensures h{t, b) < b. As we will be mainly interested 
in the distance of h to b at time t, we set 

, , ,v f b—h(t,b) if h(t, b) is well-defined 

(2-8) h(t,b):=^ Q 

This motivates to define X £ in the following way: For t £ [at, <7j_|_i), z € No, we set 
(2-9) Xt-^B^-hit-a^B^-Xl^), 

whereby X$_ := (Figure[2]). This means, for t £ [ai, ai+i) we first consider X%._. 
If B% is smaller than this value, we set X\ := B% . Else we can apply function h 
to calculate the movement of X e torwards B e . If X e reaches B e before time t, it 
remains on this level. 
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FIGURE 2. The three processes B £ (black), X e (red), and Y e 
(blue). B is displayed beneath in grey. To make the construc- 
tion clear, e is chosen large in this figure (e = 1/2). 

Lemma 2.1. Let S C [0, oo) be a compact set and e <C exp(— 7 • sup S). Then 



(2.10) 



sup 

tes 



x: 



< eKs a.s., 



whereby K,s is a finite, deterministic constant depending on S , and e' < e. 

Proof. See appendix \X\ □ 

Lemma |2~T1 shows that (X^) e>0 is a Cauchy sequence in the set of all bounded 
functions from 5 to R, equipped with the sup-norm. As this space is complete, 
(X t e ) £ >o converges. We denote the limit process by X. Equivalently, we define 

(2.11) y e (B £ (w)) := -X e (-B e (uj)) and Y(B(uj)) := -X{-B{uo)). 



3. The main theorem 

3.1. The theorem. 

Theorem 3.1. Let B, X , and Y be defined as before and let 

(3.1) 9 r :=sup{t>0: \Y t -X t \ < r} 

be the last exit time from an r-ball with respect to the \ \ ■ \ \i-norm. Then 

(1) for 7 < 1 

(3.2) (Vr > 0) 9 r = 00 a.s., 

(2) and for 7 > 1 

(3.3) (Vr > 0) 6 r < 00 a.s. 
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The theorem confirms our guess concerning 1 being a critical value for 7. For 
the critical case there is no statement at all, but as the proof of transience in the 
supercritical case seems to be sharp, our conjecture is null-recurrence if 7 — 1. 

We prove Theorem 13.11 by discretising the process Y — X. This results in a 
Markov chain, which we will examine in detail in Subsection 13.21 In 13.31 we prove 
the subcritical case by reducing it to a one-dimensional random walk problem. For 
the transient case (7 > 1) we basically use that a Markov chain is transient if we can 
find a bounded subharmonic function with respect to the generator of the chain. 
The particular theorem and its application in the proof can be found in Subsection 



3.2. Discretising the problem and facts about Markov chains. Let us look 
at the problem from another perspective. We consider the two-dimensional process 
(B e — X e , Y e — B e ), and interprete it in the following as particle moving in [0, oo) 2 . 
Observe that Y e — X e is just the sum of both coordinates. Furthermore, because 
y e — X f - can only increase at times Oi and decreases afterwards, we have 

(3-4) inf (Y<-X*) t = (Y<-X%^_. 

For all e > we define a two-dimensional Markov chain $ £ — $(i? e ) :— ($(_B e )i) igN 
with state space [0, oo) 2 , equipped with the Borel-er-algebra 5B([0, oo) 2 ), by 

(3.5) S? :={B--X%Y--B% i _, 

whereby <tq — = 0. The j-step transition probabilities from x £ [0, oo) 2 to A C 
[0, oo) 2 will be denoted by PJ(A), neglecting the index for j = 1, and the generator 
L will be given by 

(3.6) Lg(x):= [ P x (dy)g{y) - g{x) 

J[0 : oo) 2 

for suitable functions g : [0,oo) 2 — ► [0,oo). 

In the following it will be of great importance to understand, how the particle 
moves exactly, while <&\ — (x, y) jumps to (Figure [3]). At first, a jump of size 
e happens at time CTj. The position afterwards is either (x + e, (y — e) V 0) or 
{{x— e) V 0, y + e) with probability 1/2 each. Let us call this new position (x',y'). 
Before the next jump happens at time CTi+i, the particle drifts in direction of the 
origin. If it reaches one of the axis, it remains there, and only drifts forwards the 
other one, until it has reached (0,0). Thus, the coordinates of Qj+i are given by 
(/i((7j-l-i — (Tj, x 1 ), h(iji + i—<Ji,y')). Observe that <& e can only increase (in 1 1 ■ 1 1 i-sense) 
on the axes. 

Next, we need to understand the distribution of <7j-|_i — e^. Thus, we set 

(3.7) tr< := a i+1 - = inf {t > : B t = e} . 

As already suggested in the equation above, all Oi are i.i.d with support on (0, 00) 
and Ecr = e 2 . The distribution is not known explicitly, but it can be expressed as a 
series with alternating summands with decreasing absolute values (refer to section 



ESCAPING THE BROWNIAN STALKERS 



7 




I ► ^ B e - X e 

Figure 3. The particle jumps (red arrows) parallel to the level 
lines of the || • ||i-norm. In the sense of this norm it can only 
increase on the axes. The drift consists of two independent drift 
components (blue dashed arrows), orthogonal to the axes. The 
resulting drift is illustratetd by the solid blue arrow. 



C.2 in [2]). Calculating the first two summands results in 

(3.8) (l - \^ 2lt ) 

(3.9) < P(a > e) = P ( sup \B S \ < e ) 

\0<s<e / 

(3.10) < -e-* 2 /^. 

TT 

For our purposes it will be sufficient to know that both bounds are of order 
exp(-l/e). 

As we are operating on a continuous state space, the question for irreducibilty 
is a question for reaching sets instead of single states. Formally, <& e is called (p- 
irreducible if there exist a measure tp on 2}([0, oo) 2 ) s.th. 

(3.11) cp(A) > =>• P x ($ e ever reaches A) > for all x € [0, oo) 2 . 
In our case 

(3.12) P x ({0}) > for all x € [0, oo) 2 , 

because the support of a is unbounded. Thus <& e is Jo-irreducible. The existence of 
an irreducibility measure ensures that there is also a maximal irreducibility measure 
if? (compare with [7], Prop. 4.2.2) on *B([0, oo) 2 ) with the properties: 

(1) ^> is a probability measure. 

(2) <I> e is ^-irreducible. 

(3) $ e is ^'-irreducible iKVytp' (i.e. #(A) = <p'{A) = 0). 

(4) =0 => ^{{x: P x (^ e ever enters A)}) = 0. 
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(5) In our case 9 is equivalent to 

oo 

(3.13) *'(A)=J2PZ(A)2-i. 

We denote the set of measurable, ^-irreducible sets by 

(3.14) 58+ ([0, oo) 2 ) := {A e Q3([0, oo) 2 ) : 9(A) > 0}. 

Because the density of ffi+i — <ii has support on (0, oo), it is not hard to see that 

(3.15) \i(A) := Leb(A) + 6o(A) ^ => 9(A) ^ 0, 

and therefore 9 y \i, whereby Leb denotes the Lebesgue measure. 

Since 9 e is a Markov chain on the (possible) local minima of Y e — X e in the sense 
of (|3.4[) , it is obvious that transience of 9 e implies transience of Y e — X e . On the 
other hand ||$ e ||i can only increase by at most e in every step. Thus, 

(3.16) sup (y e -X e ) t <p!|| 1 + e 

t£[<T i: <T i + 1 ) 

and recurrence of $ e also implies recurrence of Y e — X e . However, observe that 
the proof of recurrence/transience for Y e — X e , e > 0, would not directly imply 
recurrence/transience for Y — X in general, because we only have convergence on 
compact sets. Thus, we will shortly argue in the end of both parts of the proof, 
why the desired result follows in our case. 

3.3. Proof of the subcritical case: 7 < 1. For the subcritical case we reduce 
the movement of $ e to a nearest neighbour random walk on the level sets 

(3.17) M(k) := {(x,y) € [0, oo) 2 | x + y = A k } , k e Z, 

of || • ||i : [0, oo) 2 — * [0,oo), and show that the probability to jump to M(k — 1) 
is larger than 1/2 + 6, 6 > 0, for small e and all k > k* for a k* e Z. Then 
it is well-known that ||$ e ||i < 4 fc infinitely often. Recurrence for 9 e follows by 
irreducibility. 

In particular, we introduce for k G Z 



(3.18) M~ (k) := { (x, y) e [0, oo) 2 | x + y < A^ 1 } , 

(3.19) M + (k) := {(x,y) e [0, oo) 2 | x + y > A k+1 } , 
and the hitting time of 9 e for a set M C [0, oo) 2 

(3.20) t € m :=mm{i : $ 2 e S M} , 



neglecting the e whenver possible. Then we have to show 

(3.21) (3k*) (Vfe > k*) lim vgPm (r e M - {k) < r^ +(fe) ) > 1/2 + 6, 6> 0. 

The proof works in four steps (Figure 2]). 

(1) We show that P m (r M -( fe ) < T M +(k)) is minimized for m* g {(4 fe , 0), (0, A k )}. 
As the model is symmetric, w.l.o.g., we may assume m* = (0,4 fc ). 

(2) We show 

(3.22) P m , (T {{x , y y. x=y} < r M+{k) ) > 1 - e^ 7 « 0.576 
as e tends to 0. 
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(3) We assume, the particle has been successful in the last step, and has reached 
(x,x) ^ M + (k). Then, in the worst case, it is at position (2-4 fe ,2-4 fe ) now, 
respectively arbitrarily close to it (as e becomes small). As the direction of 
the jumps and the drift times o~i are mutually indpendent, we can treat jump 
and drift phases independently. We will use this knowledge to determine 
the diameter of a tube around the bisector. As long as the particle is located 
within this area, it will not drift to the axis too fast. When we know the 
diameter, we can calculate the probability that the jumps do not take the 
particle out of the tube within a certain time period. Knowing this time 
and the speed torwards the origin, we can calculate, how close it gets to 
the origin before hitting the axis. 

(4) Finally we combine steps 2 and 3. It will turn out that the probability to 
stay in the tube for a certain time (step 3) can be chosen large enough such 
that it is still strictly larger than 1/2 if multiplied with the probability to 
reach the bisector (step 2). On the other hand, the time $ e stays in the 
tube, will be sufficient to reach M~(k). 

For step 1 we consider a realisation B e (uj) of the Brownian step function, and 
X e {B c (uj)) attracted to this realisation with starting distance \Bq — X$\ = d, as 
well as X e (B e (uj)), constructed like X e and attracted to the same realisation, but 
with initial distance \Bq — Xq\ = d, d > d. Because X e is Markovian, we can easily 
extend our construction of X e to initial values different from 0. Then 

(3.23) (Vt>0) (B e (uj)- X e ) f > (V(w)-X £ ) 
with equality for all t > r > 0, whereby r fulfils 

(3.24) (£»-* 6 ) =0=(b*(u>)-X< 



By symmetry the respective statement holds also for Y e ~ B e . Thus, if 
is smaller or equal in both coordinates than a copy for some time i, this 

(in)equality will remain for all times afterwards. Thus, we can conclude that for 
x < x' 

(3.25) P{xfl) ( T M-(k) < T M+(k)) > P(x',0) ( T M-(k) < T M+(k)) > 

because every realisation of B e fulfiling the event on the right side also fulfils the 
one on the left side. 

As $ e can only increase at the axes, starting it from a point inside the quadrant 
will result in a decrease of both coordinates until one of the axes is hit. But then 
(I3.25P applies, and therefore step 1 is proven. 



In step 2 we show 

(3-26) P m » (T{( x>y y. x = y } > T M +(k)) < e~ 6/7 . 

We assume m* = (0,4 fc ). The particle has two possibilities now. Either it jumps 
upwards the axis to (0, 4 fe +e), or it jumps into the quadrant to (e, 4 fc — e). Afterwards 
it drifts. In this step we will ignore the drift phase for two reasons. First, the 
change of position by jumps is of order e, while it is of order e 2 by drifting, because 
Eer = e 2 . Furthermore, the drift direction is different from the jump direction, and 
for every change of position in jump direction by drifting, there is also a drift down, 
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B e -X e 



M-(jfc) 



M(k) 



M+{k) 



Figure 4. The idea of the proof: The particle starts in M(k) 
(black line). We show that the probability to get to M~(k) (left 
dark grey area) before it gets to M + (k) (right dark grey area) is 
larger than 1/2. This probability is bounded from below by the 
product of the probability to reach the bisector (dotted line) before 
reaching M + (k), and the probability to get back to M~(k) before 
hitting the axis. In particular, we calculate the probability of a 
random walk with step size to stay in the white slot around the 
bisector. Its diameter (green line) diam^A^k+i) is a lower bound 
for the diameter of the area enclosed by g(x) and g~ x (x) (red lines). 



orthogonal to the jump direction, by the same amount at least. Thus, considering 
the drift would help us in reaching our aim to drift down. 

We introduce the following game: sitting on the axis, the particle can either 
reach the bisector or it can move up the axis by e. As the particle needs 4 fc /(2e)Q 




steps to reach the bisector, but only one step to go up, the success probability is 
small. If we should not success, we have another chance at (0,4 fe + e) (even if the 
probability for success is smaller there) and so on, until we reach M + (k). It is well 
known that the probability of an one-dimensional, symmetric random walk to reach 



Here we neglect that the expression is meaningful for integers only, because the difference will 
not play a role as e tends to zero. 
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-1 before it reaches k e N, started in 0, is given by kj(k + 1). Thus, 



(4*+ 1 -4 fe )/e-l 

(3.27) P m * {r^ x ^), x=y } > TM+[k)) = Yi 



_(A k /2 + ie)/e 
(A 

i=0 



(4 k /2 + ie)/e+l 



3-4 k /e-l 



(3.28) = n 41/2 + " 



(3.29) < 



4 k /2 + ie + e 



4 fe /2 + 3-4 fe 

(3.30) = ( I 



2e ^' 4 /e 



(3.31) -> e~ 6/7 as e tends to 0. 

For part 3 we assume that the particle has reached the bisector and is at position 
(2 • 4 k , 2 ■ 4 fc ). First, we are interested in the speed of the particle while drifting. In 
particular, we are looking for a uniform lower bound for the speed orthogonal to the 
|| • ||i-level sets on [0, oo) 2 \M + (fc). If we denote the particle's current position by 
(x, y), its speed in x-direction is given by (1 + x)~ 7 and in y-direction by (1 + y)~ 7 , 
because of equation (|2.5p . Thus, the speed orthogonal to the level sets is given by 



(3.32) «(«,„) :=V( 1+a 27 + ( 1 + 2/) 21 ■ 

Differentiation of v shows that on the set {(x,y) : x + y < 4 k+1 } the speed is 
minimized exactly on position (2 • 4 fe , 2 • 4 fe ) and amounts 

(3.33) t> mi n := V2 (1 + 2 • 4 fc )~ 7 • 

Next, let us take a closer look at the movement of the particle while drifting. 
Observe first that a drifting particle started in (x, y) will never cross the path of 
a second particle, started somewhere else, before it has hit one of the axes. This 
follows directly from our argumentation in step 1. Let us assume that x < y. By 
symmetry the other case will follow immediately. In this case, the particle will first 
hit the re-axes, and that happens at time 

(x + i)"' +1 - 1 

(3.34) t x := minit : h(x, t) = 0} = ^ , 

7+1 

which follows from the definition of h in (12. 8p . What constraints must hold for y 
such that the particle will hit the axes in M~(k)l Clearly, y must fulfil 

(3.35) h(y,t x ) < 4 fc_1 , or equivalently 



(3.36) y < ((4*- 1 + 1) 7+1 + (x+ l)^ 1 - l) 1. 

Let us denote the right side of the last inequality by g(x). By differentiation we 
immediately see that g{x) — x is a positive, strictly decreasing function, tending to 
as i becomes large. On the other hand, for starting position (x,y), y < x, the 



12 



ALEXANDER WEISS 



calculation would be the same with exchanged roles of y and x, and we would end 
up with g(y). Thus, as long as the particle starts in 

(3.37) (x,y) eA:={(x,y):(x + y< A k+1 ) A (g-\x) <y< g{x))} , 

it will first reach M~(k) and hit the axis only afterwards. This leads to the crucial 
observation: as long as the particle only jumps to positions (x, y) £ A, we do not 
have to worry that the particle will reach the axis before reaching M~{k). 
Let us define the level sets of A by 

(3.38) A t := An {(x, y) : x + y = l}. 

We can interprete A[ as a one dimensional interval or a piece of a line, and because 
g(x) — x and g~ 1 (x) — x are tending to zero, the length of this interval, denoted 
by diam(Ai), decreases as I increases. Thus, we would like to know diam(A 4 k+i), 
as it is a lower bound for all I we are interested in. Because the jump direction 
of the particle is parallel to the Ai, we can afterwards estimate, how much time 
the particle will spend in A when performing jumps. However, it is not possible to 
calculate diamlA^k+i) explicitly, but by Pythagorean Theorem, the symmetry of 
g{x) and g~ 1 (x), as well as the decrement of g(x) — x again, we have 

(3.39) diam{A Ak+ i) > V2 (g(2 ■ A k ) - 2 • 4 fe ) =: d k . 
Our ansatz is 

(3.40) d k > DA k 

for a constant D, independent of k if k is large enough. Notice that function g 
as defined in (|3.36[) is basically the || • || 7+ i-norm of (A k ^ 1 ,x) and decreases in 7. 
W.l.o.g., we may assume that 7=1. 

(3.41) V2^((4 fc - 1 + l) 2 + (2-4 fc 
easily transforms to 



, 2 



1/2 



1 -1 -1-2-4* \ > DA 



(3.42) V2D < ~ 0(4- fe ) j • 

Finally, we have to answer the question, how long do we remain in an interval of 
diameter \[2DA k , when we start in the centre and perfom a random walk with step 
size \/2e. Let us denote a standard random walk with step size 1 by R, then we are 
looking for the hitting time 

(3.43) <f (fc) := min {n : R n <£ (-DA k /e, DA k /e)} . 
It is well known that 



(3.44) Ef(fc) 



DA k 



We would like to have a lower bound for the probability that we stay in the interval 
for cE£ e (fc) steps at least, whereby c S (0,1) can be arbitrarily small. It will be 
sufficient to show that this probability tends to 1 if c goes to 0. As e tends to zero, 
Donsker's principle (see chapter 2.4.D of fi ) tells us that 

DA k - d 
(3.45) lim — — R {Di k /c) 2 t = B t , 
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whereby R is the linear interpolation of R. We define the exit time of a Brownian 
motion B from (— 1, 1) by 

(3.46) £:=inf{i:B t £(-l,l)}. 

If we use Donsker's principle we get for e tending to zero and a constant a > 

(3.47) P(C(k) <cEC(k)) = P (f < c) 

(3.48) = P (exp (-orf) > exp (-ac)) 

(3.49) < 
(3.50) 



e -ac 



cosh (v 2a) 

In line (|3.49p we have used the Markov inequality, in line (|3.50[) the explicit formula 
for the Laplace transform of £ (refer to formula 3.0.1 in [3]). As a was chosen 
arbitrary, we would like to minimize line (|3.50[) as a function of a. Differentiation 
shows that the optimizing a fulfils 

/ / — \ sinh (V2a) 
3.51) cosh ( V2a ) = \_ ' 

Using equality (|3.51| in (|3.50p results in 

(3.52) P(f (jfc) < cEf (jfc)) < -~ 



sinh (y2a) 

which tends to zero as c tends to zero. Let us call 

(3.53) p c :=P(C(k) >cW(k)) 

and observe that one can choose c s.th. p c is arbitrarily close to one. 

In step 4 we summarise the results from the steps before. When the particle starts 
in M(k), the probability to reach the bisector, before it reaches M + (k) is larger 
than 1 — exp(— 6/7) by steps 1 and 2. By step 3 we can find a c* > such that 
(1 — exp(— 6/7) )p c * > 1/2. This means, we will stay within A for c*(D4 k /e) 2 steps 
at least. As the particle drifts with a minimal speed u m i n , defined in (I3.33p . it will 
decrease its distance to the origin in terms of the || • ||i-norm by 

c*(D4 k /e) 2 

(3.54) V2(i + 2-4 fc )~ 7 ^ 

2=1 

(3.55) = c*L> 2 4 2 ' £ \/2(l + 2-4 fe )~ 7 

(3.56) = 0(V 2 -^ fc ) 

for e tending to zero. In line p.55p we have used the LLN for the i.i.d. <Ji, which 
have expectation e 2 . Thus, the distance, the particle covers, is of order ^ 2 ~i) k . 
On the other hand, the distance, the particle has to cover to get to M~(k), is by 
construction of the proof smaller or equal than 

(3.57) 4 k+1 - 4 fc - 1 = O (4 fe ) . 

Obviously (|3.56p dominates (|3.57p for 7 < 1, which finishes the proof in the sub- 
critical case for Y e — X e . 



14 



ALEXANDER WEISS 



To see that the result transfers to Y — X, we consider the process X e , constructed 
in the same way like X e but with the modified ode 

(3.58) j t f(t) = ((1 + 2c) + b - /(i))" 7 ; /(0) = 

instead of the original ode (|2.5|) . Equivalently, we define Y e (B e ) := —X e (—B e ). 
The proof easily shows that the change of the constant from 1 to 1 + 2e in (|3.58|) 
does not change the calculations or the result in an essential way (apart from 
longer equations because of the extra term 2e). Thus, Y e — X e is also recurrent 
for 7 < 1. The crucial observation is that these auxiliary processes sandwich the 
original processes: 

(3.59) Xf >X e t -t and Y t 6 ' < Y* + e 

for all e' < e. This holds due to the fact that |Xj' £ -X| e | < e and \B e - B € ' \ < e. 
Thus the difference of speed cannot be larger than It. This argument extends 
inductively to all later times <Ji. It follows 

(3.60) Y t -X t = hwy/-x t e ) 

(3.61) < Yf-XI + 2e, 
which proves recurrence for Y — X. 

3.4. Proof of the supercritical case: 7 > 1. We first define, what transience of 
Markov chains means. 

Definition 3.2. For any A C [0,oo) 2 let 

OO 

(3-62) »?A := X>{*J6A} 

i=0 

be the number of visits of <& e in A. A set A is called uniformly transient if for there 
exists M < 00 such that W,/ x ,y\(r)A) < M for all (x, y) S A. We call <I> e transient if 
there is a countable cover of [0, oo) 2 with uniformly transient sets. 

We will use the next theorem to show that <i> e is transient in the upper sense. It 
is stated as a more general result in [7j, 8.0.2(i). 

Theorem 3.3. The chain $ e is transient if and only if there exists a bounded, 
non-negative function g : [0, oo) 2 — > [0, 00) and a setB e «B+([0,oo) 2 ) such that for 
all (x,y) € [0,oo) 2 \£ ; 

(3-63) Lg(x,y)= P(x,y)(d{x,y))g(x,y) > g(x,y) 

J[0,oo] 2 

and 

(3.64) D:= J(x,y)e[0,oo) 2 



g(x,y) > sup g{x,y) } € <B+([0, oo) 2 ). 
(S,y)eB I 



Basically, we have to find a certain function g such that we jump away from the 
origin in expectation with respect to g. This must hold outside a compact set B 
containing the origin. To find a proper B we set for all z > 



(3.65) B z := G [0,oo) 2 | \\( x + l,y + 



1)IU=* 
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For g we choose 

(3.66) g(x,y):=l-\\(x+l ) y+l)\\-l 1 . 

If we can find a z, remaining finite as e tends to zero, s.th. equation (|3.63p holds 
for all (x,y) £ B z , z > z, we are done. Recall what happens in one step of $ e in 
the underlying process, described on page[H Equation (|3.63|) becomes 

(3.67) P{aedt)g{h(t,x + e),h(t,y-e)) 

i r°° 

+ 2~J P ( aedt )9{h{t,x-e),h(t,y + e)) > g{x,y) 

whereby (x,y) € B s - Because of the e-jump of B e at time a, the integral splits 
into two parts. Within both integrals the only source of randomness is a. Given 
its value, we can calculate the next position of $ e by using function h, and finally 
apply g to this value. 

Using the definition of g and observing that the integral of the density P(a £ dt) 
is one, (|3.67p easily transforms to 

1 f°° 

(3.68) -j P(aedt)\\(h(t,x + e) + l,h(t,y-e) + l)\\-l 1 
1 f°° 

+ 2J P (°£ dt )\\(Ht^-t) + hh{t,y + e)) + l\\-l x < z' 1 

As already argued, a is small, or rather we can change the upper bound of the 
integrals from oo to e at the expense of order exp(— 1/e). Furthermore, let us 
assume for the moment that x,y > 2e. As jump size and drift time are e at most 
and the drift speed is bounded from above by 1 this condition avoids that we have 
to handle cases in which the axes are reached. Observe that the only special cases 
to check later on are (0, y) and (x, 0), because we can choose for every pair x, y > 
an e > such that the condition above is fulfiled, and we let e tend to zero. Now 
we can use Taylor approximations for e and t to get 

(3.69) ~ (\\(h(t, x + e) + l, h(t, y-e) + 1)||^ 

+ \\(h(t,x-e) + l,h(t,y + e) + l)\\-l 1 ) 

(3.70) = ^((x + e + iy +1 + (y-e + iy +1 -2( 1+ l)ty^ 

+ ((x + e + 1) 7+1 + (y - e + 1) 7+1 - 2 ( 7 + 1) 

(3.71) = z- l + 2z-^+ 2 h-l((x + l)~<- 1 + (y + iy- 1 )z-^ +2 h 2 



+ (l + t)0(z-^+ 3 h 2 ). 

(3.72) ( P{a e dt)t <Ecr = e 2 

Jo 



Because 
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we can rewrite (|3.68p as 

(3.73) z- 1 + 2z^ +2 h 2 + 0(z-^ +3 h 2 ) 

< z- 1 + 1 ((x + if- 1 + (y + ly- 1 ) z-^+v 

which holds if 



> 4. 



(3.74) 7 ((* + l)^ 1 + (y + l)^ 1 ) 

Notice that equation (|3.74p is fulfiled for z large enough and 7 > 1, only. 

It remains to show the special case if x or y is zero. Because of symmetry it is suffi- 
cient to treat one of these cases. We assume x = and thus y = (z 7+1 — 1) 1 A7+ 1 )_ 1. 
Then condition (|3.67[) becomes 

-1 ^ 1 / Dfe^^llft^ ,1 , 1 Ut+ „-. ,N , 1MI-1 



(3.75) z" 1 > -/ P(aedt)\\{h(t,e) + l,h(t,y-e) + l), h _ i _ l 



P(aedt)\\(l,h(t,y + e) + l)\\;l 1 . 



Applying Taylor approximation in the same way as above results in 
(3.76) z' 1 > z- 1 ~ z-^ +2 h + 0{e 2 ) 

which is true for all 7 and arbitrary z. 

The idea, how to transfer the transient result to Y — X, is basically equal to the 
recurrent case on page 1141 This time we consider the process X e , constructed like 
X e but with the modified ode 

(3-77) j t f{t) = ((1 -2e) + b- f(t))-' ; /(0) = 

instead of (|2.5[) . Equivalently, we define Y e (B e ) := —X^(~B e ). Again, the proof is 
not essentially changed by this modifications, and thus, Y £ — X e is also transient 
for 7 > 1. Observe that the auxiliary processes are sandwiched by the original 
processes: 

(3.78) Xf <Xl + e and Y t e ' > Y t e - e 

for all e' < e. This follows from the same idea as in the recurrent case. It follows 

(3.79) Y t -X t = ]im(Y t *-Xt) 

(3.80) > Y:-Xl-2e, 
which implies the desired result. 



4. Conclusions 

In this last section we describe, what our results mean for the opinion game [3] . 
We will begin with a short description of the model. Although it is introduced in 
great generality in the original article, we will adhere to this implementation, which 
has produced interesting results in the simulations. For a deeper discussion about 
the choice of the parameters we refer to the original paper. In the second subsection 
we will point out the connections between our work and the opinion game. 
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4.1. The opinion game. Bovier et al. consider a generalised, resp. virtual, order 
book containing the opinion of each participating agent about the value of the share. 
Here the notion of value is distinguished from the one of price. While the price will 
be determined by the market and is the same for all agents, the value is driven by 
fundamental and speculative considerations, and thus, varies individually. This is a 
fundamental difference to the modelling of a classical order book. While a classical 
order book only keeps track of placed orders, the generalised order book knows the 
opinion of all market participants, independent on whether they have made them 
public. The dynamics of the model are driven by the change of agents' opinion. 

A market with N traders trading M < N stocks is considered. For simplification 
every trader can own at most one share, and furthermore, a discrete time and space 
setting is assumed. The state of trader i is given by his opinion, denoted by pi G Z, 
and the number of stocks he posseses, € {0, 1}. A trader with one share is called 
a buyer, one without a share is called a seller. The state of the order book is given 
by the states of all traders. A state is said to be stable, if the traders with the 
M highest opinions posses a share. In particular, one can fully describe the stable 
state of the order book by the price opinions p := (pi, ■ ■ ■ ,Pn) only. For stable 
states one can define an ask price as the minimum opinion of all traders possesing 
a share: 

(4.1) p a := min{p 4 : n, = 1}, 

and the bid price as the maximum opinion of all traders without a share: 

(4.2) p b := max{ K : n t = 0}. 

The current (logarithmic) price of the stock is defined by p := (p a — p b )/2. The 
update of the order book state p happens in three steps: 

(1) At time (t + 1) € No, select trader i with probability g(-;p(t),t). 

(2) The selected trader i changes his opinion to p%(t) + d, whereby d e Z has 
distribution /(•; p(t), i, t). 

(3) If p' = (pi(t), • • • ,Pi(t) + d, ■ ■ ■ ,p N (t)) is stable, then p(t + 1) = p'. Oth- 
erwise, trader i exchanges his ownership state rii(t) with the lowest asker, 
resp. highest bidder j. Afterwards, to avoid a direct re-trade, both partic- 
ipants change their opinion away from the trading price. 

The function g is defined by 

(4.3) g(i;p(t),t) := h( Pi (t) -p(t))/Z g (p(t)), 
whereby 

(4.4) h(x):=l/(l + \x\y , 7 >0, 

and Z g normalizes g, s.th. ^2^ =1 g(i' 7 p(t),t) = 1. 

The size of d is chosen from the set {—/,...,/} with probability 

(4.5) /(d;p(i),M):=7^((W)W*)) d Al) for d ? 

and f(0;p(t),i,t) = 1 — ^o<|fc|<* /(^' P(*)> *> ^ ne parameter # PijP (t) describes 
the tendency to change the opinion into direction of the price. Thus it is larger 
than 1 for Pi < p and smaller for pi > p. The second parameter, S cxt , simulates 
outer influences on the opinion change, e.g. news or rumors. This force is the same 
for all traders, but changes its strength in time. Good results were achieved by 
taking I = 4, b~ PitP (t) = exp(O.l) for buyers, and b~ Pi , p (t) = exp(— 0.1) for sellers. The 
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external influence changes its strength after independent, exponentially distributed 
times with rate 1/2000 to exp(eis' i ), whereby e$ are Bernoulli with P(ei = ±1) = 1/2 
and are Exponential with mean 0.12. Observe that in expectation the external 
force is slightly stronger than the drift to the price. 

The jump away from the trading price in the last step is implemented by setting 

(4.6) Pi (t+l)=p\t)-k, p,{t + l)=p b {t) + k 
if trader i sells a stock in this step, and 

(4.7) Pi(t + l) =p a (t) + k, Pj (t + l)=p a (t)-k 

if he buys it. In the simulations k is a uniformly distributed variable on {5, ... , 20}. 

In the simulations the price is recorded every 100 opinion updates. Thus, if 
we talk about one simulation step in the next section, we mean 100 steps of the 
underlying dynamics. 

4.2. Our result in context. Simulations show that the price process produced by 
these dynamics has some interesting properties. At first, the distribution of returns, 
that is the relative change of the price in one step, has heavy tails. Furthermore, 
the volatility, that is the average size of returns in some time interval, shows cor- 
relations on much larger time scales than the implementation would suggest. For 
the volatility of an interval of size 100, correlations after 10 4 steps can be observed. 
This is suprising, because 10 4 recorded steps are equal to 10 6 steps of the dynamics. 
But the model is Markovian and even the strength of the external influence changes 
after only 2 • 10 3 steps. 

The explanation for these observations can be found in two features of the im- 
plementation. As alreday suggested, the external force brings excitement into the 
market. Else the traders would basically perform random walks into the direction 
of the price. The returns would be much smaller, an interesting structure of the 
volatility would not exist. This coincides with the Efficient Market Hypothesis, 
because in a world without news and rumors there are no reasons for price changes. 

But the external force on its own does not explain the memory of the system in 
terms of volatility. This behaviour arises from the slower update speed of traders 
far away from the current price. This mechanism makes sure that the system 
remembers price changes on large time scales. If we observe an order book state in 
which a group of traders has a large distance to the current price, we can conclude, 
the price must have been in the region of the traders before, as it is very unlikely 
that a whole group of traders has moved against its drift. Furthermore, after fast 
price movements the distance between ask and bid price, called gap, is larger than 
average and needs some time to recover. In these periods the market is illiquid and 
a small number of trades can move the price a lot, which results in an increased 
volatility. Increased volatility after large price movements is a well observed feature 
of real world markets. 

Thus the connection of update speed and distance to the price is of paramount 
importance for the model. Indeed, the larger 7 is chosen in formula (|4.4p the better 
the just explained phenomena can be observed. However, a larger 7 contains the 
risk of instablity of the whole system. It turns out that once the gap has exceeded a 
certain size (depending on 7), it cannot recover anymore and the two groups, buyers 
and sellers, drift away from each other. Then the price waves between these groups, 
driven by two traders, one from each group, which were able to get away and now 
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Figure 5. Screenshots of the virtual order books after 428500 sim- 
ulation steps for 7 = 1.5 (left) and 7 = 1.6 (right) with same initial 
conditions and same realisation of external influences. Observe the 
different distances between buyers (green) and sellers (red) and the 
different behaviour of the price processes (blue box). 



1 

- 








y 

1/ / 


/ 

/ ; 







III 

1 



25B+Q5 5e+05 7.5e+05 le+06 1.25e+06 1 5e+06 



lc+(J5 2c+U5 3c+U5 4c+U5 5e+05 



Figure 6. The left graph shows the gap of the system for different 
7. While it is stable for 7 = 1.5 (black lower graph), it increases for 
7 = 1.6 (red) and 7 = 1.7 (grey). However, if the system is started 
with 7 = 1.5 but with an artificially enlarged gap, it also increases 
(black increasing graph). The convergence to a value below 2000 is 
due to a restriction of the state space in the numerical simulations. 
The right graph shows the stable resp. unstable behaviour for 
7 = 1.5 (black) and 7 = 1.6 (red) in terms of the price process. 



basically move according to the external drift without any resistance by surrounding 
traders. For 7 > 1.6 this happens quite fast while the model has remained stable in 
simulations over several days for 7 = 1.5 (Figure [5]). On the other hand, if we start 
a simulation already with a large gap and 7 = 1.5, also this system is not able to 
recover. As a large gap size will eventually reached by randomness, it is justified 
to talk about a metastable behaviour. In Figure [6] we illustrate these statements 
with a sample. Instead of recording the difference between ask and bid price, we 
have taken the distance between the 950th and the 1050th trader ordered by their 
opinions (i.e. the buyer with the 50th highest opinion and the seller with the 50th 
lowest one), because traders close to the price suffer much more fluctuations than 
agents with some distance. In this sense our choice represents the majority of the 
traders. 

In the situation when the trader groups have already a large distance from each 
other, the two traders in between, and also the price, perform basically a random 
walk. Especially, when the two traders are close to the middle in between both 
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groups, their probability to move is almost 1. In this case our model with a Brow- 
nian motion as driving force offers a reasonable approximation for the behaviour of 
the system. Thus, our results give few hope that any simulation with 7 > 1 will be 
stable forever. But for 7 < 1 the memory effect producing all the statistical facts is 
too small. However, as already mentioned, the model seems to be stable on a large 
time scale for 7 = 1.5. This and also the sharp threshold between 1.5 and 1.6 are 
not understood. More research is neccessary here. 

Besides these findings the three particle model introduced in this paper has its 
qualities on its own. As a simple model for longtcrm investors, this easy setting 
already exhibits an interesting and non-trivial longterm behaviour. As a logical 
next step it will be interesting to see, how the results change if we substitute the 
Brownian motion by a Levy process, which is much more realistic for price process 
on stock markets. 



Appendix A. Proof of Lemma 12711 
We turn to the proof of Lemma |2~T1 

Let S C [0, 00) be a compact set and e -C exp(— 7 • sup S). Then 



(A.1) 



sup 



< eKs a.s., 



whereby K$ is a finite, deterministic constant depending on S, and 



Because S is compact, w.l.o.g., we may assume S = [0,t*] for some < t* < 00. 
Remember that the jump times of B e were denoted by a e in (|2.1|) . and the time 
between two jumps by <r e in (|3.7|) . Furthermore, 



(A.2) 



B e -B e 



We denote the distance of X e to B t by 
(A.3) d t := B%. 

and the distance to X e by 
(A.4) At := XI 



< e. 



XL 



XI 



always meaning er with respect to e. We would like to maximize A 2 , thus, we 
assume that B e has jumped upwards at o\. Then d\ = e and |Aj| < e. We first 
assume that Ax is positive. By definition of A and of h in (|2.6p . 



(A.5) 



A, 



(X, 



,T, X aiJ 



x k-xh 



X', 



x% 



IA~2l 



(A.6) < h(a 1 ,d 1 )-h(a 1 ,d 1 +A 1 +e) + A 1 

(A.7) IP h(a 1 ,d 1 + A 1 +e)-h(a 1 ,dx)-e. 

Remember that h is basically defined as 

, . . - , . , V(7+l) 

(A.; 



MM) = ((d + i) 7+1 -(7 + i)*)" 



- 1. 
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As the distance will not increase anymore, once X e has hit B e , we get an upper 
bound for a\\ 

<di + 1 % 1 7+1 - 1 

(A.9) h{a u d x ) > tr x < \ +l • 

Because d\ — e we have a\ < e. As di, Ai, e, and o\ are small in comparison to 1, 
we apply Taylor twice to line (|A.7[) and get 

(A.IO) A 2 < (l-( 7 + l)a 1 )^ /(7+1) (A 1 +e)-e 

(A.ll) = Ai + 7 (Ai + e) <7i 

(A.12) = A 1 (l + 7 e). 

With the same argumentation we can conclude that 

(A.13) A l+1 < Ai(l + 7e), 

and thus, 

(A. 14) X\,-Xi = At. /e 

(A.15) < A 1 (l+ 7 e) t * /£ 

(A.16) -» ee 7t *. 

On the other hand, if JT e > AT e , basically the same idea applies: the distance grows 
the quickest, if one of the processes always stays close to its attracting process s.th. 
it has drift speed 1. Now, if X e increases with speed 1 (as a worst case assumption), 
a c = e and we end up with the same calculation as before. 

It should be mentioned that our estimations are rough, as we do not consider 
the structure of Brownian paths, but only the worst case of all continous paths. 
However, uniform convergence on compact intervals is the best one can get and 
every improvement would only change the constant K$. 
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